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Scanning tunneling images of carbon nanotubes frequently show electron distributions which break 
the local sixfold symmetry of the graphene sheet. We present a theory of these images which relates 
these anisotropics to the off diagonal correlations in the single particle density matrix, and allows 
one to extract these correlations from the observed images. The theory is applied to images of 
the low energy states reflected at the end of a tube or by point defects, and to states propagating 
on defect free semiconducting tubes. The latter exhibit a novel switching of the anisotropy in the 
tunneling image with the sign of the tunneling bias. 



PACS: 72.80.Rj, 73.40.Gk, 61.16.Ch 

Scanning tunneling microscopy and spectroscopy is a 
powerful tool for studying the structural and electronic 
properties of carbon nanotubes at the atomic scale. Sev- 
eral experimental groups have reported tunneling images 
of isolated single wall carbon nanotubes [0,|| and of tubes 
packed into bundles or "ropes" 0. In some cases these 
measurements have allowed a direct determination of the 
diameters and wrapping vectors for the tubes and these 
observations, combined with scanning tunneling spec- 
troscopy, have confirmed the idea that the semiconduct- 
ing or conducting behavior of a tube is controlled by its 
wrapping vector [^J|]. 

However STM images of these systems obtained at low 
bias voltages contain a number of surprising features. At 
low bias voltages the images rarely display the full sixfold 
lattice symmetry even when the underlying graphene lat- 
tice is undistorted. Instead, these images frequently con- 
tain a broken symmetry in the form of "striped" patterns 
in which maxima in the electron density are observed in 
bond chains which spiral around the tube &-|3] . In some 
images superlattice structures are present with a period 
commensurate with but larger than that of the underly- 
ing graphene sheet [Q. Moreover, energy resolved images 
of short tubes show standing waves characteristic of in- 
dividual eigenstates which also have a period longer than 
that of the graphene lattice g . 

In special cases broken translational or rotational sym- 
metries obtained in an STM image can be attributed to 
asymmetries in the tunneling tip g]. In this Letter we 
point out that asymmetric images are expected even in an 
ideal tunneling experiment and contain important infor- 
mation about the low lying electronic states in these sys- 
tems. The asymmetries are interference patterns which 
are sensitive to the coherence between the "forward" and 
"backward" moving electronic states propagating on the 
tube walls. This can arise from backscattering from tube 
ends, from various defects on the tube walls, and even 
from propagation in a translationally invariant potential 
on a semiconducting tube. We show that these inter- 
ference patterns in the tunneling images are fingerprints 
which directly probe the off diagonal correlations in the 
single particle density matrix. We provide a theory for 



extracting these correlations from the observed images. 
This effect is illustrated with several examples of the tun- 
neling densities of states calculated for defect free tubes 
and for tubes with point defects. 

The low lying electronic states on a carbon nanotube 
are derived from the propagating states near the "Fermi 
points" of an isolated graphene sheet located at the 
Brilloun zone corners shown in the lower panel of Fig. 
1 [05§]- There are two inequivalent points which we 
will refer to as K = Kq and K'=-i^o- Wrapping the 
graphene sheet into a cylinder requires that the elec- 
tronic wavefunctions satisfy periodic boundary condi- 
tions ^(r + Tm,n) = ^(0 where Tm,n = mri + nT2 gives 
the wrapping vector around the tube circumference ex- 
pressed in terms of the two primitive graphite translation 
vectors n = (1/2,^3/2) and r2 = (-1/2,^3/2). When 
mod(7Ti — n, 3) = the Bloch waves exactly at the K and 
K' points are allowed quantized waves on the tube. This 
leads to the metallic band structure in Fig. 1(a). In con- 
trast, when mod(m — n, 3) = (1, 2) the allowed quantized 
waves do not intersect K and K', which leads to a semi- 
conducting gap in the electronic spectrum as shown in 
Fig. 1(b). 

The spectrum of Fig. 1(a) describes the propagating 
modes of a defect free conducting tube. However, these 
waves can be reflected from tube ends or from defects 
along the tube. The interference between the forward and 
backward moving waves produces a spatial modulation 
of the charge density. Since the carbon nanotubes have 
two forward and backward moving channels (associated 
with the K and K' points shown in Fig. 1 the resulting 
interference patterns have a particularly rich structure. 
The coherent superposition of the forward and backward 
moving components of the scattering states produces off 
diagonal correlations in the density matrix at energy E^ 



pME) = {^PiSiE - H)#). 



(1) 



where H is the Hamiltonian and a is a four component 
index specifying the left and right moving bands at the 
K and K' points. 

To derive the local tunneling density of states from the 
density matrix in equation (1) we represent the Bloch 



waves V^aC^) as a sum atomic orbitals centered on sites 
Tm in cells T^ 

Mr) = ^7mc.e^^^-^-/(r - r^ - f,). (2) 



where ^ma are the amplitudes for the Bloch state on 
sites m. These Bloch waves can be represented as an 
expansion in reciprocal lattice vectors 



^a = J2^rnae-'^^'+'^"'>-^-F{\k^+Gn\)e' 



(fca+G„)-r 
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For tunneling from a tip which is smooth on a scale of 
the atomic spacing, F{q) decreases rapidly for large q. 
In the following we will truncate this expansion, keeping 
only G's in the lowest "star" of ka + G. This becomes 
exact when the STM tip is sufficiently high above the 
surface [B. Including higher Fourier components does 
not significantly change our conclusions. We also assume 
the tip is isotropic, so within the first star F{\k -\- G\) is 
independent of G. 

The local density of states at energy E can be ex- 
pressed 



Kr,E) = C(^>a/3(^)#(r-). 



(4) 



It is useful to characterize the tunneling image in terms 
of its longest wavelength Fourier components. Coupling 
between bands at the same K point leads to images with 
the periodicity of the lattice. These are described by 
Fourier components in the first star of reciprocal lattice 
vectors qu = (5^, indicated in Fig. 1. On the other 
hand coupling between the two K points leads to modu- 
lations with a V^ X V^ superlattice, which are described 
by the "a/3" star, q/^. = Ki. The discussion is further 
simplified by projecting onto the "triangular harmonics" 
defined in each star to be 



n=0 



/"dVe-*«''"V(f) 



(5) 



with p = 1 or a/3 and m = —1,0,1. Combining (1- 
5) yields a simple expression for the Fourier components 
ppm{E) in terms of the density matrix papiE). 

We first consider the effects of reflection either from 
the end of the tube or from an impurity. The reflection 
of waves associated with the K and K' points is charac- 
terized by a 2 X 2 matrix of complex reflection ampli- 
tudes. This matrix contains three independent ampli- 
tudes labeled r^, r^ and r^ in Fig. 1(a). These describe 
respectively large momentum scattering from K ^ K' ^ 
K' ^ K and small momentum scattering at the K and 
K' points. The equality between the two amplitudes 
described by r^ follows from time reversal symmetry. 
These reflection amplitudes depend on the detailed struc- 
ture of the scatterer, although for special high symme- 
try scatterers its form can be constrained by symmetry. 
In the following, we consider an infinitely long tube, so 



that the density matrix is diagonal in the basis of scat- 
tering states, which are a superposition of incoming and 
reflected waves. 

The large momentum scattering amplitudes r^ and r& 
illustrated in Fig. 1 produce a modulation of the TDOS 
in the ^/?> star. For a conducting tube we find that the 
TDOS measured at energy E and at a distance x from a 
point scatterer has the Fourier coefficients [|lO[ 



p^o(E) = \N{E){ne 



2iQx _ ^*^ — 2iQx 
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PV3±ii^) = |Af(£;)e±*^(rbe2*'5" + r*e-2'<3^) 
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where N{E) is the density of states, Q = E/fivp^ vp is 
the Fermi velocity and 6 is the chiral angle which ori- 
ents the zigzag bond direction with respect to the tube 
axis (i.e. ^ = defines a tube with an armchair wrap- 
ping). The dependence of the TDOS on r^ and r^ is 
most clearly seen for tunneling at low bias {E ^ {)). The 
Q dependence in equation (5) arises from the fact that 
the scattered states at nonzero energy are not located 
precisely at the Fermi points. 

As an example in Fig. 2 we display the calculated 
TDOS at E=0 for the values r^ = r^ = 1 on an armchair 
tube. The scattering amplitudes produce a a/3 x a/3 mod- 
ulation of the tunneling image which in which the bond 
charges are enhanced in a superlattice of bonds oriented 
along the circumferential direction of the tube. Similar 
a/3 X a/3 modulations occur in the presence of impurities 
on the surface of graphite |ll|. Those patterns follow 
from the two dimensional scattering between the K and 
K' points in the graphite plane. 

The small momentum backscattering amplitudes Vm 
produce cell periodic modulations of the TDOS which can 
nonetheless break the rotational symmetry of the image. 
These effects are produced by a modulation of the Fourier 
components of the TDOS in the first star of reciprocal 
lattice vectors Gn shown in Fig. 1. Using the expansion 
in triangular harmonics in equation (5) we find [ffol 






(7) 



Interesting structure in the TDOS is produced by the 
imaginary part of pi±i. This generally occurs for any chi- 
ral tube with 6> 7^ 0, but it can also occur for a nonchiral 
armchair tube, ^ = 0, when the q ^ backscattering am- 
plitude rm develops a nonvanishing imaginary part. This 
leads to the interference pattern shown in Fig. 2(b), in 
which a bond density wave is deflected into a spiral pat- 
tern around the axis of the armchair tube. We find that a 
reflection amplitude with this symmetry can be produced 
by any point defect or end cap which breaks the two sub- 
lattice symmetry of the underlying graphene sheet. 

Semiconducting tubes with mod(7Ti — n, 3) 7^ have a 
gap in the low energy spectrum as shown in the middle 
panel of Fig. 1. This gap arises from a coherent super- 
position of forward and backward moving components 
which is required for the wavefunction to satisfy periodic 



boundary conditions around the tube waist. Exactly at 
the band edges one obtains a perfect standing wave which 
contains an equal admixture of forward and backward 
propagating waves. Note however that the backscatter- 
ing responsible for these states does not result from re- 
flection from an isolated point defect, but instead arises 
from the presence of a "mass" operator in the low energy 
Hamiltonian H which preserves the lattice translational 
symmetry, and breaks its rotational symmetry. Remark- 
ably, the symmetry of this "mass" term depends sensi- 
tively on the wrapping vector and allows one to distin- 
guish mod(7Ti — n, 3) = 1 from mod(7Ti — n, 3) = 2 tubes. 
To do this we deflne the chiral index s = (_i)mod(m-n,3) 
Then we flnd that the tunneling density of states for a 
semiconducting tube with chiral index 5, chiral angle 6 
and gap 2 A is given by jl^] 






(8) 



In Fig. 3 we display tunneling densities of states cal- 
culated for the band edge states {E = ±A) on tubes 
with wrapping indices [m,n] = [12,7], [14,6] and [17,0] 
(zigzag). These three tubes have a common chiral index 
5=1 and chiral angles which vary between 8.6° and 30°. 
The band edge states imaged at positive and negative 
bias have complementary structures, that is the super- 
position of the two images gives an image with perfect 
sixfold symmetry although each image separately breaks 
the sixfold symmetry. For chiral angles near 0^0 the 
tunneling images consist of a series of complementary spi- 
ral stripes. However the TDOS changes smoothly as a 
function of chiral angle, so that as one approaches the 
zigzag structure the negative energy states are enhanced 
in a pattern of isolated bonds in the structure, while the 
density for the positive energy states is conflned to a 
connected zigzag bond chain. Note that the symmetry 
breaking terms in pi±i depend on the product of the chi- 
ral index and the energy. Thus for a [12,8] tube which 
has a chiral index -1 the symmetries of the positive and 
negative energy solutions are reversed. 

Fig. 3 shows that for semiconducting tubes the tunnel- 
ing images obtained near the band edges break the local 
point symmetry of the graphene sheet, and that for a 
given tube the sign of the symmetry breaking is switched 
by reversing the bias of the tunneling tip. Observation 
of such a reversal in the tunneling image would provide a 
striking identiflcation of the chiral index s of a semicon- 
ducting tube, even when the wrapping indices [n,m] can 
not be resolved. Moreover this reversal would clearly dis- 
tinguish this effect from a tip artifact ^ . We should note 
that several experiments have already suggested that a 
tube in contact with a conducting substrate can be doped 
"p" type [Q and so one needs to correct such a measure- 
ment for the offset in the chemical potential for any such 
unintentionally "doped" tube. One also needs to ensure 
that the tunneling is carried out in an energy window 
where only one azimuthal subband is accessible, since the 
symmetries of the tube eigenstates alternate in successive 



conduction or valence subbands, tending to suppress the 
anisotropy in the tunneling images. 

The data presented in Figures 2 and 3 demonstrate 
that scanning tunneling images of carbon nanotubes are 
very sensitive to the nodal structure of the underlying 
electronic states [^. Indeed at low energy these images 
are probing the internal structure of individual (or at best 
a small number) of electronic eigenstates on the tube sur- 
face. The appearance of these broken symmetry patterns 
in a tunneling image does not imply large structural per- 
turbations to the covalent graphene network. Indeed the 
data presented here are calculated for structures which 
are all unstrained and perfectly locally sixfold symmet- 
ric. These density patterns do require coherence between 
forward and backward propagating waves. In fact a quan- 
titative analysis of these images can be used to extract 
off diagonal correlations in the density matrix in equation 
(1) responsible for these patterns. These data can then 
be used to extract the scattering matrix, which charac- 
terizes the internal structure and symmetry of various 
scattering centers, and to identify the wrapping vector 
for semiconducting tubes. Finally, we remark that the off 
diagonal correlations responsible for these patterns could 
also arise due to interactions with a substrate, interac- 
tions between tubes in a rope or due to electron-electron 
interactions. 
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FIG. 1. Low energy electronic spectra for a conducting 
tube (top) and for a semiconducting tube (middle). The lower 
panel shows the momenta in the first two stars in reciprocal 
space which describe the local tunneling density of states in 
an STM experiment. 





FIG. 3. Complementary tunneling densities of states for 
the two band edges of three semiconducting chiral tubes with 
the same chiral index. The left panel gives the image for tun- 
neling into a band edge state at negative energy (E = —A), 
and the right panel gives the corresponding image at positive 
energy (E = A). The images are for tubes with wrapping 
indices [12,7] (top), [14,6] (middle) and [17,0] (bottom). In 
each case the image is rotated so that the tube axis is di- 
rected along the horizontal. 



FIG. 2. Local tunneling densities of states for a [10,10] tube 
in the presence of reflection from point defects. The left panel 
shows a tunneling image which is modulated in a Vs x a/3 
pattern and in the right panel shows a "primitive" 1x1 spi- 
ral striped pattern which breaks a mirror symmetry of the 
armchair tube. 



